Constraining the Symmetry Parameters of the Nuclear Interaction 
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We show combining nuclear mass data with experimental information from neutron skins, heavy 
ion collisions, giant dipole resonances and dipole polarizabilities, and theoretical calculations of 
neutron matter, lead to stringent constraints on the symmetry properties of the nuclear interaction 
near the nuclear saturation density. Furthermore, these constraints are remarkably consistent with 
inferences from astrophysical observations of neutron stars. The concordance of experimental, the- 
oretical and observational analyses suggests that neutron star radii, in the mass range 1 Mq - 2 
Mq, lie in the narrow window 11 km < i? < 12 km. 

PACS numbers: 26.60.Kp, 97.60. Jd, 21.65. Cd 
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Introduction.- The nuclear symmetry energy, and in 
particular its density dependence, has been the focus 
of much recent activity due to its importance to astro- 
physics ^. The expected neutino signal from gravita- 
tional collapse supernovae within our Galaxy is sensitive 
to the symmetry energy both before core bounce [2] and 
at later times from proto-neutron stars j3]. The radii 
of neutron stars are closely correlated to the pressure of 
neutron star matter in the range Us — Sn^ @] , and there- 
fore to the density derivative of the symmetry energy 
at the nuclear saturation density n^. Neutron star radii 
strongly affect both moments of inertia [5] and the ex- 
pected gravitational signals from mergers involving neu- 
tron stars 6 . The composition of matter at densities 
above Us , and the existence of neutrino processes which 
can rapidly cool neutron stars, depends on the density 
dependence of the symmetry energy [7] , as do the prop- 
erties of neutron-rich nuclei and reaction rates involved in 
the astrophysical r-process. The symmetry energy con- 
trols properties of neutron star matter below Ug , such as 
the core-crust transition density in neutron stars which is 
important in pulsar glitches |5], crust relaxation in cool- 
ing [S] and accreting [T^ neutron stars, and astroseismol- 
ogy from giant magnetar flares [TT]. 

Nuclear Symmetry Energy. - The nuclear symmetry en- 
ergy is alternatively expressed in one of two ways: as 
the difference S between the energies of pure neutron 
and isospin symmetric matter, S{u) = e(u, 0) — e(u, 1/2), 
where e(u, x) is the energy per baryon of nuclear matter 
with baryon density n — uug and proton fraction x, or as 
the coefficient 5*2 of a quadratic expansion of symmetric 
matter in the neutron excess. 



e[u, x) = e{u, 1/2) + S-2.{u){\ - 2xf + . . 



(1) 



It is often assumed that 5' ~ S'2 but this need not be true. 
While the kinetic contributions from non-interacting nu- 
cleons lead to a small difference, 

8Sk^n{u) = Sk^n{u) - S2.k^n(u) ~ 0.72^2/3 McV, (2) 

there is no guarantee that potential energy contributions 
are similarly small. Many Skyrme interactions were con- 
structed with the property that their potential energy 



contributions to bS^ot vanish, but some other Skyrme in- 
teractions have unphysical neutron matter properties at 
high densities due to large potential energy contributions. 
It is customary to define the liquid droplet symmetry 
energy parameters S^ , L and Ksy„i : 
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While these parameters are not useful in describing the 
nuclear symmetry energy far from the saturation density, 
most astrophysical applications are most sensitive to the 
symmetry energy near the saturation density. There are 
limited constraints on Kgy^ from measurements of the 
giant isoscalar monopole resonance [T^], isotopic trans- 
port ratios in medium-energy heavy ion collisons |13) , and 
neutron skin data from anti-protonic atoms |14| , but ex- 
traction of the symmetry contribution to the total incom- 
pressibility of nuclei and its separation into bulk {Kgy^) 
and surface contributions is still problematic. 

In contrast, there are several experiments leading to 
effective constraints on S^ and L, including nuclear mass 
fitting, neutron skins, dipole polarizabilities, dynamics in 
heavy ion collisions, giant and pygmy dipole resonances, 
and isobaric analog states. 

Nuclear Mass Fitting.- The most accurate and least 
ambiguous data are nuclear masses. However, optimizing 
the values of symmetry parameters from experimentally- 
measured masses cannot uniquely constrain them due 
to their significant correlation [15]. This can be eas- 
ily observed by referring to the standard liquid drop 
model [T5I in which the symmetry energy of an iso- 
lated nucleus is Edrop.A = I^{SvA — S'^A^/'^) where 
I ^ {N - Z)/{N -(- A) is the neutron excess. The bulk 
symmetry parameter S^ is identical to that appearing in 
Eq. (|3|. 5s is the surface symmetry energy parameter 
which depends on the behavior of S within the surface 



and thus S'„, L and Kg 



The correlation follows from 



minimizing the difference between the model and exper- 
imentally measured symmetry energies, i.e., minimizing 



x'^E^iEe 



Edrop.iY /{J^f^"^), where N is the total 



number of nuclei and a a nominal error. A contour of x^ 
one unit above the minimum value represents a Icr confi- 
dence ellipse. The properties of the ellipse are determined 
by the second derivatives of x^ ^t the minimum. 
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where Xvs = d^X^ /dS^dSs, etc. In this simple example, 
these derivatives depend only on Ai and /i, not the loca- 
tion of the x^ minimum. The confidence ellipse in Ss — 
Sy space has orientation a — (l/2)tan~^ \'^Xvs/{Xvv — 
Xss)\ — 9.8° with respect to the Ss axis, with error widths 
<^v = ^/{x~^)vv — 2.3cr and a^ = \/ix^'^)ss ^ 13.2ct 
where (x~^) is the matrix inverse, and a correlation co- 
efficient r = Xvs/^/XvvXss ^ 0.997. 

It is necessary to transform the correlation between Ss 
and Sy into one between L and Sy. Treating the nuclear 
surface as a plane-parallel interface, to lowest order Ss 
can be expressed as the ratio of two integrals over the 
density profile u = n/us of a symmetric system [1] : 






Es /o u^'^jSylS^ju) - l]f{u)-^'^du 
2 J^u^f^f{uy/^du 



(5) 



where f{u) = e(u, 1/2) + B, B = -e(l,l/2) ~ 16 
MeV is the bulk symmetric binding energy at n^, and 
Es ~ 19 MeV is the liquid drop surface energy param- 
eter (also estimated from fitting nuclear masses). The 
transformation depends on 5*2 and e and is generally com- 
plex, but an analytical result is found if we approximate 
S2{u) ~ Sy+L{u-l)/3 and f{u) ~ i^o(M- 1)^/18, where 
Ko — 240 MeV is the incompressibility: 

Ss/Sy ~ 135^, [1 - /3tan-i r'] /C^K), (6) 



where (3 = y^3Sy/L - 
finds Ss/Sy 



1. 



With L/Sy = 1, 3/2, 2, one 
0.7, 1.2, 1.8 which shows the basic trend. 
A much better approximation comes from semi-infinite 
Thomas- Fermi calculations [17]: 
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Transforming to Sy — L space, one finds a confidence 
ellipse with Uy ~ 2.3cr and ol — 23tT. 

In practice, the liquid droplet model [12], which ac- 
counts for varying neutron/proton ratios within the nu- 
cleus, represents an improved treatment of the total nu- 
clear energy. Its symmetry energy is £^ld,a — /^>5'„^(1-|- 
SsA~^/^ / Sy)~^ and therefore predicts a linear correla- 
tion between Ss/Sy and Sy. The same methodology as 
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FIG. 1: Comparison of confidence ellipses for nuclear mass 
fitting. The solid ellipse is the EDF \a result of Kortelainen 
et al. [20] and the dashed ellipse is the liquid droplet model 
result assuming a — 0.5 MeV. The quantities Svo and Lq 
(plotted as circles) are the minimum x^ values and (e^)^" 
are the RMS errors of the respective fits. 



for the drop model can be used to determine the confi- 
dence ellipse. One ultimately finds ay ~ 2.7 MeV and 
ctl — 37 MeV if ct = 0.5 MeV (these values now de- 
pend on the location of the minimum, Syo ~ 29.6 MeV 
and Lq ~ 47 MeV, with shell and pairing effects taken 
into account). As shown in Fig. [ll the centroid and 
orientation are in remarkable agreement with a recent 
Hartree-Fock fit to experimental masses and radii us- 
ing a parametrized energy-density functional (EDF) [20] , 
which finds Syo :^ 30.5 MeV and Lq ~ 45 MeV. The 
RMS error of the EDF fit is approximately 1.45 MeV, 
with ay = 3.1 MeV and a^ — 40 MeV; these values are 
similar to those of the liquid droplet fit if ct « 0.5 MeV. A 
parametrized finite-range Thomas-Fermi (FRTF) model 
for nuclei produces a nearly identical result [IT . The cor- 
relation of Sy and L is approximately 97% in the EDF fit, 
smaller than the correlation of the liquid droplet fit since 
the EDF model was also constrained by nuclear radii. 

The high degree of correlation of the symmetry param- 
eters established from nuclear mass fitting has already 
been used to constrain equation of state tables [TH] , but 
can be combined with data from other experiments to 
produce even greater constraints. 

Other Experimental Constraints.- The liquid droplet 
model implies that other nuclear observables will be re- 
lated to the surface and volume symmetry coefficients. 
An example is the neutron skin thickness of neutron- 
rich nuclei. Neglecting Coulomb effects, the difference 
between the mean neutron and proton surfaces in the 
liquid droplet model 19 is 
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which is primarily a function of Sg/Sy. It should then 
form a nearly orthogonal constraint to that from nuclear 
masses. Values of tnp have been measured, typically with 
30-50% errors, for a few dozen nuclei. A recent study of 
Sn isotopes [H], where differential measurements with 
fixed Z can reduce errors, found the correlation between 
Sv and L plotted in Fig. [2J It is nearly orthogonal to the 
mass-fit correlation (the 95% confidence ellipse from Ref. 
PP] is shown in Fig. [2]). First results [H] from the PREX 
experiment to measure the neutron skin of ^"^Pb are too 
imprecise to be definitive as they imply 35 < L/MeV < 
262 for reasonable values of Sy. The central value of L is 
discrepant with other results for neutron skins, but the 
uncertainties are still very large. 



droplet parameters. The centroid energy is 
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where k is an enhancement factor arising from the ve- 
locity dependence of the interaction and < r^ > is the 
mean-square charge radius. The factor m < r^ > /{1 + k) 
does not significantly vary among interactions for a given 
nucleus. As a result, a series of measurements of the 
GDR among various nuclei should result in a correlation 
between Ss/Sy and Sy that is similar to but slightly less 
steep than that from nuclear masses. The liquid droplet 
model predicts the orientations of the correlations in the 
Ss/Sy — Sy plane will have the ratio 
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FIG. 2: Summary of constraints on symmetry energy param- 
eters. The filled ellipsoid indicate joint Sv — L constraints 
from nuclear masses [2D]. Filled bands show constraints from 
neutron skin thicknesses of Sn isotopes [21], the dipole po- 
larizability of Pb [25], giant dipole resonances [171 126] . and 
isotope diffusion in heavy ion collisions [57]. The hatched 
rectangle shows constraints from fitting astrophysical M — R 
observations |37) . The two closed regions show neutron mat- 
ter constraints (H is Ref. [33] and G is Ref. [3^). The white 
area is the experimentally-allowed overlap region. 



The hydrodynamical model of Ref. [53] showed that 
both the dipole polarizability and the giant dipole res- 
onance centroid energy are closely connected to liquid 



where <> refers to an average over the measured nu- 
clei. Using a variety of Skyrme functions, Trippa, Colo 
& Vigezzi [26] show the GDR in 2°*Pb imphes that S2{u) 
evaluated at the density n = 0.1 fm^'^ is limited to 23.3 

< 'S'2.o.i/MeV < 24.9. While the connection between 
'S'2,0.1 and Sy and L is model dependent, bounds for a 
range of parametrizations [17] are indicated in Fig. [2] 
The correlation band has a shallower slope, compared to 
that from masses, as expected from the droplet model. 

Tamil et al. [M] established the limit 23 < i/MeV < 
54 from measurements of the dipole polarizability of 
^•^^Pb. A systematic analysis incorporating model de- 
pendencies [25] revised this range to 37 < L/MeV < 60, 
and is shown in Fig. [2J 

Constraints from heavy ion collisions include those 
from isospin diffusion [27] . shown in Fig. pi and from 
multifragmentation in intermediate-energy heavy ion col- 
lisions [32]. The latter implies 40 < L/MeV < 125 but is 
not displayed in Fig. [2] for clarity. 

Not included in Fig. [2] are results of a study of pygmy 
dipole resonances [55] in ^^Ni and ^'^^Sn which imply 31 

< Sy/UeV < 33.6 and 49.1 < L/MeV < 80. Daoutidis & 
Goriely [29j argue that theoretical and experimental un- 
certainties prevent them from being effective constraints. 
The results of Reinhard & Nazarewicz [30| confirm pygmy 
resonance energies lack significant correlation with L. 

Danielewicz & Lee [JT] proposed isobaric analog states 
as constraints on symmetry parameters, leading to Sy ~ 
33 MeV and S^ ~ 96 MeV, which is equivalent to L ~ 113 
MeV. However, the data suggest there is a significant 
curvature contribution; re-fitting with curvature yields 
Sy ~ 30 MeV and S^ ~ 45 MeV, equivalent to L ~ 50 
MeV, but with considerable errors [T7] . A more detailed 
analysis is necessary to produce better constraints. 

In summary, the various experimental results are con- 
sistent with each other, and they jointly confine the 
allowable ranges of symmetry parameters to the small 
white overlap region in Fig. [2] 



Neutron Matter and Astrophysical Studies.- Recent 
studies of pure neutron matter have employed chiral la- 
grangian [33] and quantum Monte Carlo [34] techniques 
using realistic two- and three-nucleon interactions. The 
energies and pressures of neutron matter at Ug , if quartic 
and higher order symmetry energy contributions can be 
neglected, are e(l, 0) = Sy — B and p(l, 0) — Lns/3. De- 
spite the different theoretical approaches, and consider- 
ing physically- motivated ranges for interaction strengths, 
the neutron matter predictions for the symmetry parame- 
ters are consistent with each other and with experimental 
studies, as shown in Fig. [2] 

Neutron star radii are highly correlated with neu- 
tron matter pressures near rig [1]. Quiescent neutron 
stars in globular clusters |35| provide information con- 
cerning their apparent angular radius, Roo/D, where 
Roo = i?(l - 2G'Af/i?c2)-i/2 and M and R are the 
stellar mass and radius. Although the distances D are 
only approximately known, these sources are believed to 
have well-understood H-dominated, low magnetic field, 
atmospheres. Type I X-ray bursts on neutron star sur- 
faces [55] provide, in addition to angular sizes, measured 
maximum fluxes thought to be near the Eddington limit 
FEdd — cGM/{kD'^) (neglecting redshift) at which radi- 
ation pressure balances gravity {k is the opacity). For 
these sources, major uncertainties include distances as 
well as atmospheric properties. Despite estimated errors 
in M and R of order 20% or more for each source, a 
Bayesian analysis [37] of an ensemble of both types of 
sources implies 11 < R/km. < 12 for masses between 1.2 
and 2 M©, and 40 < L/MeV < 60 (see Fig.[2|. Interest- 
ingly, pulse-shape modeling [38] of accreting millisecond 
pulsars produces a virtually identical constraint. 

Conclusions.- There is emerging concordance among 
experimental, theoretical and observational studies for 
values of the nuclear symmetry energy parameters and 
neutron star radii. This will reduce uncertainties from as- 
trophysical simulations of supernovae, neutron star merg- 
ers, proto-neutron stars and cooling neutron stars. 
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